Complex singularities of the critical potential in the large- limit 
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We show with two numerical examples that the conventional expansion in powers of the field for 
the critical potential of 3-dimensional 0{N) models in the large- limit, does not converge for values 
of (/>^ larger than some critical value. This can be explained by the existence of conjugated branch 
points in the complex (j)^ plane. Pade approximants [L + 3/L] for the critical potential apparently 
converge at large 4>^ . This allows high-precision calculation of the fixed point in a more suitable set 
of coordinates. We argue that the singularities are generic and not an artifact of the large- A'^ limit. 
We show that ignoring these singularities may lead to inaccurate approximations. 



I. INTRODUCTION 

Since the early days of the renormahzation group (RG) 
method [1], 3-diniensional scalar models have been iden- 
tified as an important laboratory to discuss the existence 
of non-trivial fixed points and the large cut-off (or small 
lattice spacing) limit of field theory models. In the case 
of iV-components models with an 0{N) invariant La- 
grangian, the RG transformation becomes particularly 
simple in the large- limit [2]. The construction of the 
effective potential for these models is discussed in Refs. 
[3-5]. Later, motivated by perturbative results indicating 
the existence of an UV stable tricritical fixed point for 
N large enough [6], a new mechanism allowing sponta- 
neous breakdown of scale invariance and dynamical mass 
generation was found in the large- iV limit [7] . In the fol- 
lowing, we call this mechanism the "BMB mechanism" . 
It was argued [8] that the BMB mechanism is compatible 
with a zero vacuum energy and a better understanding 
of this question might suggest a solution to the cosmo- 
logical constant problem. Spontaneous breaking of scale 
invariance is also discussed [9] with related methods in 
four-dimensional models of clear interest in the context 
of particle physics. However, doubts were cast [10] about 
the fact that the BMB mechanism is generic and it is 
commonly believed that it disappears at finite N. 

In this article we report results which force us to re- 
consider the way we think about non-trivial fixed points. 
We usually think of the RG flows as taking place in a 
space of bare couplings or more generally in a space of 
functions. The necessity for this more general point of 
view appears quite clearly in exact renormahzation group 
equations [11]. Unfortunately, it seems impossible to de- 
cide a-priori which space of functions should be consid- 
ered to study the RG flows. It is clear from perturbation 
theory that near the Gaussian fixed point, low dimen- 
sional polynomial approximations of the local potential 
should be adequate. However, it is not clear that this 



kind of approximation should be valid far away from the 
Gaussian fixed point and in particular near non-trivial 
fixed points. 

In the following, we concentrate on the non-trivial fixed 
point found numerically in the case iV = 1 by K. Wil- 
son [1]. This fixed point is located on a hypersurface of 
second order phase transition which separates the sym- 
metric phase from the broken symmetry phase. In the 
following we call this fixed point the Heisenberg fixed 
point (HFP for short) as in Ref. [10]. It should not be 
confused with the fixed point relevant for the BMB mech- 
anism and which is not studied in detail here. The main 
result of the article is that the bare potential correspond- 
ing to the HFP has singularities in the complex plane and 
that ignoring these singularities may lead to inaccurate 
approximations. These claims are based on explicit cal- 
culations performed in the large- limit for two 0{N) 
invariant models reviewed in section II. These two mod- 
els are: 1) a model with a fc^ kinetic term together with 
a sharp cut-off, the sharp cut-off model (SCM) for short; 
2) Dyson's hierarchical model (HM) [12,13]. 

Before entering into technical details, three points 
should be clear. First, all the results presented here are 
based on the analysis of long numerical series and no at- 
tempt is made to give rigorous proofs. Second, in order 
to understand some of the statements made below, the 
reader should be aware that even though, at leading order 
in the large- A'^ approximation, the critical exponents take 
A'^-independent values, the same approximation provides 
finite A^ approximate HFP which are A'^-dependent. A 
more precise formulation of this statement can be found 
in Sections II and VI. Third, we only work in 3 dimen- 
sions. The precise meaning of this statement for the hi- 
erarchical model is explained at the end of section II. 

In section HI, we review the basic equations [2,10] 
defining the HFP for the SCM. We then show that the 
definition can be extended naturally for the HM. The cor- 
rectness of this definition is verified later in the paper. In 
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section IV, wc present the methods used to calculate the 
critical potential expanded as a Taylor series in (jp' . The 
coefficients of this expansion are called the critical cou- 
plings. The main conclusion that we can infer from our 
numerical results is that the Taylor series is inadequate 
for large values of (jP. First of all, one half of the critical 
couplings are negative. If we truncate the Taylor series 
at an order such that the coefficient of the highest order 
is negative, we obtain an ill-defined functional integral. 
In addition, the absolute value of the critical couplings 
grows exponentially with the order and the expansion 
has a finite radius of convergence. Consequently, the idea 
that the critical potential associated with the HFP can 
be approximated by polynomials should be reconsidered. 

It is nevertheless possible to define the critical theory 
by using Fade approximants for the critical potential. In 
section V, we show that sequences of approximants con- 
verge toward the expected function in a way very remi- 
niscent of the case of the anharmonic oscillator were the 
convergence can be proven rigorously [14]. In addition, 
the zeroes and the poles of these approximants are lo- 
cated far away from the real positive axis and follow pat- 
terns that strongly suggest the existence of two complex 
conjugated branch points. 

The complex singularities of the critical potential 
should not be interpreted as a failure of the RG approach 
but rather as an artifact of the system of coordinates 
used. In section VI, we present consistent arguments 
showing that in a different system of coordinates [15,16], 
the function associated with the HFP is free of singu- 
larities. In this system of coordinates, finite dimensional 
truncation is a meaningful procedure which, in the case of 
the HM, allows comparison with independent numerical 
calculations at finite N . An example of such a calculation 
is presented in the case N = 5. 

In section VII, we discuss the errors associated with 
two approximate procedures that can be used to deal 
with the singularities. The first procedure which is jus- 
tified in the context of perturbation theory and does not 
require an Tinderstanding of the singularities, consists in 
truncating the potential at order (0^)''. The second pro- 
cedure consists in restricting the range of integration of 
(jy^ to the radius of convergence of the critical potential. 
If the range of integration is large enough, this second 
procedure generates small errors [17]. As far as the calcu- 
lation of the HFP in the system of coordinates of Section 
VI is concerned, both procedures have a low accuracy for 
both models considered. In the conclusions, we explain 
why we believe that the singularities persist at finite N 
and we discuss implications of the existence of these sin- 
gularities for other problems. 



II. MODELS 

We consider lattice models defined by the partition 
function 



Z{J) = ll d^<l>,e-'+^^'-'<^- , (1) 

with 

S=-\Y.k^.y$y + Y.Vo{(t>l) . (2) 

xy X 

We use the notation = (l>x-<Px and A^y is a symmetric 
matrix with negative eigenvalues, such as discrete ver- 
sions of the Laplacian. For the simplicity of the presen- 
tation, we will assume that Axy ~ 0. If it is not the 
case, one can always subtract the zero mode from A and 
compensate it with a new term in V. 
Defining the rescaled potential 

Vo{X) = NUoi^) , (3) 

and performing a Legendre transform from the source J 
to the external classical field one can show that in 
the large N Hmit [10] that = 2dVeff/d<pl obeys the 
self-consistent equation 

2Ul,{cf>l + U{M')) = M\ (4) 

where JaIM^) is the one-loop integral corresponding to 

the quadratic form A and a mass term M^. The prime 
denotes the derivative with respect to the 0{N) invari- 
ant argument. The explicit form of / for the two models 
discussed in the following are given in Eqs. (6-7). Pre- 
cise definitions of 0^ and the effective potential Vef / are 
given in [10] . 

Up to now, all the quantities introduced are dimension- 
less. They can be interpreted as dimensionful quantities 
expressed in cut-off units. Let us consider two models, 
the first one with a rescaled potential Uq, a UV cutoff 
A and a quadratic form A and a second model with a 
rescaled potential Uo,s. a UV cutoff A/S and a quadratic 
form As- For £> = 3 and in the large- iV limit, the 
two models have the same dimensionful zero-momentum 
Green's functions provided that: 

U'o,si^') = (5) 

In two special cases, the dimensionless expression for the 
one-loop diagram is independent of the cut-off. In other 
words, /a = Jas = f and the fixed point equation be- 
comes very simple [2,10]. 

We now discuss the two models where this simplifica- 
tion occurs. In the SCM, A becomes k'^ in the momentum 
representation (Fourier modes) . The momentum cutoff is 
sharp: < 1 (in cut-off units). This is why we call this 
model the sharp cutoff model. The non-renormalization 
of the kinetic term is justified in the large- limit [10]. 
For this model, 



2 



f d^k 1 

J\k\<i (27r)^ + z 



(6) 



By construction [13], the kinetic term of Dyson's hier- 
archical model (HM) is not renormalized and we have 



fHM{z) = 



°° o-n-1 



^ b(c/2)" + z 



(7) 



with c = 2^~^/^ and h=^^. The inverse temperature /3 
and the parameter c appear in the hamiltonian of the HM 
in a way that is explained in section II of Ref. [18]. The 
parameter c is related to the dimension D by considering 
the scaling of a massless Gaussian field. In the following 
wc will consider the case c = 2^/'' (D = 3) exclusively. In 
addition, /3 will be set to 1 as in other fixed point calcu- 
lations [18]. Different values of /? can be introduced by 
a trivial rcscaling. Note also that the cutoff camiot be 
changed continuously for the HM, because the invariance 
of / is only valid when we integrate the degrees of free- 
dom of a the largest momentum shell (corresponding to 
the hierarchically nested blocks in configuration space) 
all at the same time. For the HM, the density of sites is 
reduced by a factor 2 at each RG transformation. The 
linear dimension ( "lattice spacing" ) is thus increased by 
by a factor 2^/^ and the cutoff decreased by the same 
factor. Consequently, for the HM, Eq. (5) should be 
understood only with S = 2'/^ = 2'/^ and q integer. 



III. THE HFP 

In this section, we review the construction of the HFP 
for the SCM, and we show that the construction can be 
extended in a natural (but non-trivial) way for the HM. 
The first step consists in finding all the fixed points of 
the RG Eq. (5). Following Refs. [2,10], we introduce the 
inverse function: 



(8) 



and the function H(z) = F{z) — f{z) , where the one- 
loop function / has been defined in the previous section 
for the two models considered. With these notations, the 
fixed point equation corresponding to Eq. (5) is simply 



H{z) = SH{z/S^) . 



(9) 



For the SCM, S is allowed to vary continuously in Eq. 
(9) and the general solution is 



F{z) = fscMiz) + Kz'/^ . 



(10) 



For the HM, S can only be an integer power of 2^/^ and 
the general solution has an infinite number of free pa- 
rameters: 



F{z)=fHM{z)+J2^'i^'^'^'"' ' 



(11) 



with 



37r 
h^2 



13.6 , 



(12) 



and q runs over positive and negative integers. The only 
restriction on the constants K and Kg is that F should 
have a well defined inverse which is real when F{= 0^) 
is real and positive. 

It is clear from Eqs. (6) and (7) that for both mod- 
els f{z) has singularities along the negative real axis and 
that, in general, F{z) cannot be defined for z real and 
negative. This imposes restrictions on the choice of the 
constants K and Kg. For instance, in the case of the 
SCM, when K takes a large positive value, it is impos- 
sible to reach small values of F = (f)^ when z > and 
the fixed point has no obvious physical interpretation. 
However, there is a special positive value of K for which 
the singularity of fscM is exactly canceled and an ana- 
lytic continuation for z < is possible. Its exact value 
can be calculated [10] by decomposing fscM into a regu- 
lar part fscM,reg. and a singular part fscM,sing.- Using 
elementary trigonometric identities, one finds 

fscM,reg.{z) = ^(1 + z^/^Avcta.n{z^/^)) , (13) 



and 



fsCM,sing.iz) 



.1/2 



47r 



(14) 



F reduces to 



Consequently, if we choose K = 

fsCM,reg.- 

A more detailed analysis [10], shows that this value of 
K is the only positive value of K for which U' can be 
defined for any real positive (fp. On the other hand, for 
negative K, one obtains a line of fixed points ending (for 
K = 0) a,t the fixed point relevant for the BMB mecha- 
nism. Given the isolation of the fixed point with K = j^, 
it is easy to identify it with the HFP. We denote the cor- 
responding inverse function by -FscmC-^) = fsCM,reg.i.z). 
As promised this function is analytical in a neighborhood 
of the origin and has a Taylor expansion: 



Fscm{z) = ^{1 + z- 



z^ z^ 



(15) 



This expansion has a radius of convergence equal to 1 
due to a logarithmic singularity at z = —1. However, 
as wc will see in section IV, this expansion allows us to 
construct an inverse power series and Uq. 

In the case of the HM, the decomposition into a reg- 
ular and singular part is more tedious. Fortunately, this 
problem is a particular case of a problem solved in sec- 
tion 5 of Ref. [19] where Eq. (5.6) with A = (?,B = c'^ 
and f[z) = G{z/b) yields 



fHM,sing. (^) = ~ ^ X] 



(z/b) 



1/2+iqui 



sin(7r(l/2 -|- iquj)) 



(16) 
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with b and c defined in seetion II. 

If we compare this expression with the general solution 
of the fixed point equation (11), we see that in both ex- 
pressions, the power 2;^/^+*'''^ appears for all positive and 
negative integer values of q. There exists a unique choice 
of the Kg in Eq. (11) which cancels exactly the singu- 
lar part of /hm- We call the corresponding fixed point 
the HFP of the HM. The numerical closeness with the 
finite N HFP discussed in section VI confirms the valid- 
ity of this analogic definition. We call the corresponding 
inverse function F'^j^. Using Eq. (5.5) of Ref. [19], we 
find 

1 °° / —z \ ' 1 

Fhuiz) = fHM,reg. = ^J2(~r) I _ c'^l-l ' (^'^) 



bt 




This expansion has a radius of convergence bc^ = 
2.7024 ... for the choice of parameters used here. 

It is possible to check the accuracy of the expansion 
given in Eq. (17) by using the identity F^f^[z) = 
fHM{z) — fHM,sing.- Notc that the two terms of the 
r.h.s. cannot be defined separately on the negative real 
axis. On the real positive axis, fHM,sing. is dominated 
by the q = term. Numerically, 

Stt 

^° = 46^ = '-^^°^^^---- (^^^ 

The terms with g = ±1 produce log-periodic oscillations 
of amplitude 1.7 x 10""'^^. The terms with larger j^l have 
a much smaller amplitude. These findings are consis- 
tent with the log-periodic oscillations found numerically 
in HT expansions [20,19]. The oscillatory terms are very 
small along the positive real axis. However, in the com- 
plex plane, if we write z = rc*^, the amplitude is multi- 
plied by e~*'^^ which compensates the suppression of the 
denominators in Eq. (16), if — > -|-7r(— tt) when q < 
{q > 0). In conclusion, along the real positive axis, we 
can use the approximation 

F^M{z)^fHM{z) + Koz'/\ (19) 

with an accuracy of 18 significant digits, but this approx- 
imation is certainly not valid near the negative real axis. 



IV. CALCULATION OF THE CRITICAL 
POTENTIAL US 

In the previous section, we have provided power series 
for the inverse function F(z) corresponding to the HFP 
of the SCM and the HM. We can use these series to de- 
fine F{z) on the negative real axis. In both cases, as we 
move toward more negative values oi z, F becomes zero 
within the radius of convergence of the expansion. The 
situation is illustrated in Fig. 1 for the HM. 



FIG. 1. Fff]^{z) versus z 

Numerically, we have i^^^^(— 1.5107 .. .) = and 
FgQj^{— 0.694:8 .. .) = 0. We then reexpand the series 
about that value of z (which corresponds to F = (p"^ =0) 
and invert it. The resulting scries is an expansion of 2Uq^ 
in (jy^. After integration, and up to an arbitrary constant 
uo, we obtain a Taylor series for the critical potential Uq 
corresponding to the HFP. We denote the expansion as 

00 

W) = E""('^')"- (20) 

n=0 

The precise determination of the zero of F is obtained 
by Newton's method with a large order polynomial ex- 
pansion. This expansion is then rcexpandcd about the 
zero and the large order coefficients in the original ex- 
pansion have an effect on the low order coefficients of the 
reexpanded series. We have checked that the order was 
sufficiently large to stabilize the results presented here- 
after. 
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10 20 30 40 50 

Order n 



FIG. 2. Natural Logarithm of the absolute value of the co- 
efficients Un of the critical potential Uq defined in Eq. (20) 
for the SCM (filled squares) and the HM (empty circles). 
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The absolute values of the first 50 coefficients of both 
models are shown in Fig. 2. In both cases, it appears 
clearly that the absolute value grows at an exponential 
rate. Linear fits of the right part of Fig. 2 suggest a 
radius of convergence of order 0.11 for the SCM and 2.5 
for the HM. The signs of both series follow the periodic 

pattern: H h H h H for the SCM and + H for 

the HM. This suggests singularities in the complex plane 
at an angle with respect to the positive real axis (fc = 
1, 3, 7, 9) for the SCM and along the imaginary axis for 
the HM. This analysis is confirmed by an analysis of the 
poles of Pade approximants presented in the next section. 



V. PADE APPROXIMANTS OF Uq 



At this point, our series expansion of the critical po- 
tential does not allows us to define the critical theory as a 
functional integral. As cj)^ exceeds the critical values esti- 
mated in the previous section, the power series is unable 
to reproduce the expected function Uq. The situation is 
illustrated in Fig. 3 for the HM. 



u 



Parametric 
Pade [4/1] , [5/2] 

Series (50) 




FIG. 3. C/^(0^) for the HM with a parametric plot (filled 
squares), the series truncated at order 50 (thick solid line) and 
Pade approximants [4/1] (thin line slightly above the squares) 
and [5/2] (thin line closer to the squares). The constant has 
been fixed in such a way that the value at the minimum is 
zero. 

The numerical values of Uq in Fig. 3 have been cal- 
culated using a parametric representation (with z as the 
parameter). We have calculated pairs of values 



(21) corresponds to a pair ((^^, [/q ((/)^)) with the arbi- 
trary constant in Uq hxed in such a way that Uq van- 
ishes at its minimum. We have calculated F* by us- 
ing the independent but approximate Eq. (19). As ex- 
plained in section III, the approximate expression is only 
valid for z real and positive and should give 18 correct 
significant digits. In Fig. 3, we have used the values 



z = 2m = 0, 0.25,0.5, 



This is why the filled squares 



only appear when the derivative of Uq is positive. Unlike 
Eq. (17) which has a radius of convergence 1, the approx- 
imate expression Eq. (19) remains valid for large positive 
values of z. It is thus possible to check if Pade approxi- 
mants can be used to represent the critical potential be- 
yond the radius of convergence of its Taylor expansion. 
Fig. 3 shows that low order approximants are close to 
the parametric curve. As the order increases, the curves 
coalesce with the parametric curve and a more refined 
description is necessary. 

In Fig. 4, we give the accuracy reached by various ap- 
proximants for the HM with a broad range of 0^ (more 
than 4 times the radius of convergence). As the order of 
the approximants increase the accuracy increases but at 
a rate which is slower for larger values of 0^ . The figure is 
very similar to sequences of Pade approximants obtained 
for the ground state of the anharmonic oscillator (see Fig. 
1 of Ref. [17]), where the convergence can be proven rig- 
orously [14]. Note that the slow convergence at large (jP 
is not a serious problem, since the contributions for large 
(jy^ are exponentially suppressed in the functional inte- 
gral. The choice of [L + 3/L] approximants is discussed 
in more detail below. Up to now, we only discussed the 
HM. Following the same procedure for the SCM, we ob- 
tain very similar figures (with a diff'erent scale) which 
we have not displayed. 
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for various real positive values of z. A simple graphi- 
cal analysis performed by representing [/g as a surface 
on Fig. 1 shows that each pair of values in in Eq. 



FIG. 4. Number of correct significant digits obtained with 
Pade approximants [9/6], [13/10], [17/14] and [21/18] for var- 
ious values of 0^ for the HM. 

The singularities oi Uq in the complex 0^ plane can 
be inferred from the location of the zeroes and poles of 
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the Pade approximants. As L becomes large, regular 
patterns appear. Examples are shown in Fig. 5 for the 
SCM and Fig. 6 for the HM. In both cases, the zeroes 
and poles approximately alternate along two lines end- 
ing where singularities were expected from the analysis 
of coefficients in Section IV. This pattern suggests [21] 
the existence of two complex conjugated branch points 
at the end of these lines. 
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FIG. 5. Real and imaginary parts of the roots of the 
denominator (filled squares) and numerator (crosses) of a 
[26/23] Padc approximant for the SCM The solid circle has 
a radius 0.11 and the two solid lines make angles ±^ with 



respect to the positive real axis. 
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FIG. 6. Real and imaginary parts of the roots of the de- 
nominator (filled circles) and numerator (crosses) of a [26/23] 
Fade approximant for the HM. The solid circle has a radius 
2.5. Two roots farther away on the imaginary axis and one 
root farther away near the negative real Eixis are not displayed. 



The choice of [L + 3/L] approximants is easily justi- 
fied for the SCM. At large |2;|, fscM{z) oc 1/z and and 



F^CAii^l^ 47^'^'- For large \cj)% Sn^^^ and 

Uq ~ ^(0^)'^. Consequently a [L + 3/L] approximant 
should have the correct asymptotic behavior. More pre- 
cisely, if ai_|_3 and are the leading coefficients of the 
numerator and denominator of a Fade [L + 3/L] respec- 
tively, we expect that when L is large 



bL 



3 



Defining a quantity 



Et 



3aL 



+3 



(22) 



(23) 



that measures the departure from the expected asymp- 
totic behavior, we see from Fig. 7 that as L increases, 
the discrepancy diminishes exponentially. 



O 

5 




FIG. 7. Log\EL\ versus L. 

In the case of the HM, the situation is more intricate. 
From Eq. (19), wc may be tempted to conclude that the 
two cases are similar. Unfortunately, Eq. (19) is a real 
equation not a complex one. In the complex plane, the 
terms with q ^ become important near the negative 
real axis and no simple simple limit as in Eq. (22) ap- 
plies. However, if we need Uq only along the real positive 
axis. Fig. 4 justifies the use of the [L + 3/L] sequence of 
approximants. 



VI. THE HFP IN A CONVENIENT SET OF 
COORDINATES 

As explained in the introduction we can think that the 
RG flows move in a space of functions. The system of 
coordinates for this space can be chosen in a way which 
is convenient to make approximations. A particularly 
convenient system of coordinates for the HM consists in 
considering the Fourier transform of the local measure of 
integration [15,16]. In this system of coordinates and at 



6 



leading order in the 1/N expansion, the HFP for a given 

iV reads: 

R*{k) OC j d^^e-i'^'-NU^oi'^VN)+ik.$ _ (24) 

The quadratic term proportional to 6 is due to the fact 
that the quadratic form A for the HM has a zero mode. 
We then Taylor expand 



i?*(fc) = l + 5^a„(fc2)" 



(25) 



though iV = 5 is not a large number, the first coefficients 
obtained in the leading order in the approximation 
coincide with about two significant digits with the accu- 
rate values found numerically with N = 5. As the order 
increases, the accuracy degrades slowly. This is explained 
in more detail below. However, the ratios of successive 
coefficients still follows closely the asymptotic prediction 
obtained from Eq. (27). This strongly suggests that the 
{(p"^)^ asymptotic behavior of the critical potential per- 
sists at finite A''. 



and consider the the our new set of coordinates. 

The advantage of this representation is that it is possible 
to make very accurate calculations by using polynomial 
approximations [15,16,18] of the infinite sum in Eq. (25). 
In this section and the next section, we discuss the de- 
tails of the calculations for the HM. The case of the SCM 
shares many similarities with the HM and is discussed 
briefly at the end of each section. 

We have performed a numerical calculation of the a„ 
of the HM using Eq. (24) in the particular case = 5. 
The study of the ratios of successive coefficients displayed 
in Fig. 8 indicates that the |a„| decay faster than 1/n! 
and that R*{k) is analytical over the entire complex 
plane in contrast to C/g (</>^) which has a finite radius of 
convergence in the complex (fy^ plane. 

The good convergence of R*(k) can be explained by 
an approximate calculation. The </) integral that is per- 
formed in the calculation of the a„ has a positive inte- 
grand with a peak moving to larger values of |(^| when 
n increases. For sufficiently large values of n, we can re- 
place Uq by its asymptotic behavior on the positive real 
axis which can be derived from the approximate Eq. 
(19) for the HM: 



1 r 



R*{k) ~ J d^c, 



-(i/(m^Kl))(4>^Y+ik.,t> 



(26) 



With this approximation, the a„ can be expressed exactly 
in terms of gamma functions and a simple calculation 
yields 



, (6A^^A:g)i/3r((jV + 2n)/6) 

_i ~ An{n - 1 + N/2)T{{N + 2(n - l))/6) 



(27) 



Note that there are no free parameters in this formula. 
Fig. 8 shows that Eq. (27) is a very good approximation 
of the ratios obtained numerically from Eq. (24). 

We have also calculated the a„ corresponding to the 
HFP for A" = 5 using the numerical method developed 
in the case A^ = 1 in Ref. [18] and which can be extended 
easily for arbitrary A. In brief, it consists of finding the 
stable manifold by fine tuning the temperature and then 
iterating the RG transformation in order to get rid of 
the irrelevant directions. This procedure is very accu- 
rate and completely independent of the approximations 
made in this article. Remarkably, we found that even 
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FIG. 8. Ratios of successive coefficients for tfie HM, using 
the leading order Eq. (24) (stars), the asymptotic formula 
Eq. (27) (continuous line) and the numerical fixed point 
(empty circles). Same results for the SCM: leading order 
(filled squares) and asymptotic (dashed line). In all cases, 
A = 5. 

Except for the comparison with independent numeri- 
cal calculations at finite A^, the same calculations can be 
performed for the SCM with minor changes (6^0 and 
Kq ^ K). The results are also shown in Fig. 8 where one 
can see that the agreement with the asymptotic formula 
is very good even at low order. 



VII. DISCUSSION OF ALTERNATE 
PROCEDURES 

In section V, we have shown that the Fade approxi- 
mants provide accurate values of far beyond its radius 
of convergence. In order to estimate the error on the new 
coordinates a„ due to the use of approximants for Uq, we 
can vary the range of integration and change the approx- 
imants. For instance the values of o„ of the HM used in 
Fig. 8 have been calculated using a range of integration 
1^1 < 20 and a [26/23] Fade approximant. For the values 
of n considered here, changing the range of integration 
has effects smaller than the errors due to numerical in- 
tegration (which has an accuracy of about 11 significant 
digits in our calculation) provided that we include values 
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up to \(t>\ ~ 4.9. Restricting the range of integration to 
smaller values produces sizable effects. As an example, 
the small effects due a restriction to \(p\ < 4.4 are shown 
in Fig. 9. Similarly, the values of a„ are not very sensi- 
tive to small changes in the Fade approximants. Sizable 
effects are obtained by changing the order of the numera- 
tor and denominator by approximately 10. For instance, 
the effects of using a [14/11] approximant are shown in 
Fig. 9. 
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FIG. 9. Number of significant digits common with our best 
estimate for the a-n obtained for the HM from Eq. (24) with 
n (the order) going from 1 to 20. The alternate procedure 
are the truncation at order (</>^)^ (filled circles), the N = 5 
accurate numerical result (filled square), no Fade approxi- 
mants but a truncation of the range of integration close to 
the radius of convergence (empty circles), a restriction of the 
range of integration for 4> < 4.4 (stars), and a [14/11] Fade 
(diamonds) . 

Having demonstrated that we can calculate the first 
20 coefficients a„, at leading order in the 1/A'' expansion, 
with at least 10 significant digits, we can now discuss the 
errors associated with other procedures mentioned in the 
introduction. The first procedure consists in tnmcating 
Uq keeping only the terms up to order (0^)"^. This is pro- 
cedure inspired by perturbation theory amounts to keep 
only the relevant and marginal directions near the Gaus- 
sian fixed point. From Fig. 9, we sec that this procedure 
generates errors which are of the same order as the errors 
due to the use of the leading 1/N approximation. Con- 
sequently, this procedure is quite unsuitable to study the 
correction to this approximation. Slightly better results 
are obtained by keeping as many terms as possible in the 
expansion (up to 50 in our calculation) but restricting the 
range of integration in such way that we stay within the 
radius of convergence. Given the rescaling of Eq. (3) this 
means that for = 5, we need to restrict the integration 
to |(/)| < \/5 2.5 ~ 3.54 which is substantially smaller 
than the acceptable field cutoff 4.9 mentioned above. As 
one can see from Fig. 9, this creates errors which are be- 



tween one and two orders of magnitude smaller than the 
1/A^ corrections. This is better but it compares poorly 
with what can reached with Fade approximants. 

Again, except for the comparison with independent nu- 
merical calculations at finite N, the same calculations can 
be performed for the SCM with minor changes. Results 
very similar to those shown in Fig. 9 for the HM can be 
produced. Since it contains essentially the same infor- 
mation, it has not been displayed. It should however be 
noted that the number of significant digits obtained with 
the two alternate procedures are lower than in the case 
of the HM. In the case of the truncation of the range of 
integration, we need to restrict to < \/5 x 0.11 ~ 0.74 
while a range of about 2 is required in order to obtain an 
accuracy consistent with the method of numerical inte- 
gration. 



VIII. CONCLUSIONS 



We have shown in two differents models where the crit- 
ical potential can be calculated at leading order in the 
1/A^ expansion that these potentials have finite radii of 
convergence due to singularities in the complex plane. 
Do such a results persist at finite A''? In the case of the 
HM, the behavior of the ratios at finite N shown in Fig. 
8 strongly suggests that at large real positive c6^, the 
critical potential still grows like (0^)^. Can an infinite 
sum converging over the entire complex plane have this 
kind of behavior? This is certainly not impossible (e.g., 
{(j)"^)^ + e""^ ), however it requires cancellations that we 
judge unlikely to happen. Consequently, we conjecture 
that the singularities observed are generic rather than 
being an artifact of the large- A" limit. 

We have observed that in a system of coordinates 
where the HFP can be approximated by polynomials, 
the procedure which consists in considering bare poten- 
tial truncated at order ((/>^)^ describes the HFP with a 
low accuracy. We are planning to investigate if similar 
problems appear near tricritical fixed points. In particu- 
lar, reconsidering the RG flows in a larger space of bare 
parameters may affect the generic dimension of the in- 
tersections of hypersurface of various codimensions and 
help us finding a more general realization of spontaneous 
breaking of scale invariance with a dynamical generation 
of mass. 

Our results have qualitative similarities common with 
those of Refs. [22]: we found some "pathologies" which 
force us to look at the RG transformations in a more 
open-minded way. We are planning [23] to compare in 
more detail, the leading order results presented here with 
finite A^ results, as suggested in Ref. [24] for the local po- 
tential approximation. Another issue regarding the 0{N) 
models and which would deserve a more detailed inves- 
tigation is the question of first order phase transitions 
[25,26]. 
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